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Abstract 

We explore the possibility to trap polynitrogen clusters inside Cqo fullerene cage, opening a 
new direction of developing nitrogen-rich high energy materials. We found that a maximum of 13 
nitrogen atoms can be encapsulated in a Cqo cage. The nitrogen clusters in confinement exhibit 
unique stable structures in polymeric form which possess a large component of (~ 70-80%) single 
bond character. The iV n @C6o molecules retain their structure at 300K for n<12. The Mulliken 
charge analysis shows very small charge transfer in N@Cqq, consistent with the quartet spin state 
of N. However, for 2<n<10, charge transfer take place from cage surface to N n compounds and 
inverse polarization thereafter. These nitrogen clusters when allowed to relax to N2 molecules 
which are triply bonded are capable of releasing a large amount of energy. 
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1. INTRODUCTION 



In recent years, interest in polynitrogen or polymeric nitrogen has drawn considerable 
attention not only because of its theoretical significance but also its possible application 
as high energy density material (HEDM) 

HQ. 

Polynitrogen molecules are formed by a 
combination of lower order bonds and while decaying into dinitrogen (A^) molecules, they 
release enormous amount of energy and are environmentally safe. 

Over last few decades, there have been consistent efforts to predict new exotic forms 
of all-nitrogen molecules using various experimental and theoretical techniques 0, 0, 5), Q]. 
The higher members of nitrogen cluster family (N n for n>3) are unstable in the free space 
and therefore present a challenge to their synthesis. In 1999, Christe et al.\]\ synthesized 
successfully a N§ salt N + AsF~ which is third readily accessible homonuclear polynitrogen 
species. Subsequently, the samegroup reported the stability of N^SbF^ up to 70° C and its 
relative insensitivity to impact [8]. This discovery opened a new dimension to explore neutral 
polynitrogen compounds and basis for the first synthesis of stable nitrogen allotropes. 

Nitrogen rich compounds such as nitrogen hydrides (N n H m ) (n>m) could also be used 
as starting materials to produce N n species, since N 4 H + } N 4 H 6 and N 4 H 4 0,Q,0 species 
are stable and their subsequent deprotonation may lead to the formation of neutral N4. 
Similarly, the oxidation of nitrogen oxide N 4 which is stable and has been observed ex- 
perimentally [jjj], would also lead to the formation of N 4 . Further, it has been suggested 
that higher N n could be synthesized via direct or indirect excitation of electronic states in 
liquid and solid nitrogen followed by collision addition of a ground state N and excited N n _i 
or N n _ 2 thereby forming higher N n compounds. For example, by combining different N n 
moieties such as N*. and A3 might lead to formation of N* species Q3I and NrN^ may 

nn n 

form iV 10 |14j . I15I ]. In 2002, Cacace et al. [16[ demonstrated the existence of the tetran- 
itrogen molecule (N 4 ) as a metastable species whose life time exceeds 1 /^s at 298K. The 
identification of N4 represents the first addition in nearly half a century to the family of 
the polynitrogen molecules. Hammerl and Klapotke Yj\ also reported a combined theoret- 
ical and experimental NMR study on nitrogen rich compounds. At the most fundamental 
level, a molecule will exist only if it has high enough energy barrier that keeps it away from 
dissociation. 

This possibility has been explored theoretically. Various novel polynitrogen species have 
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been studied in variety of forms such as cyclic, acyclic or caged conformations. The species 



investigated include ionic clusters 



dies 



19 



20| | . cylinders [21], cages 



26( | and helices 27] . Isomers of smaller systems N 7 



22 



23. 



24 



251 ]. nanonee- 



28|, A^ 10 [29j and N 12 [30j have 



also been reported theoretically in free space. Majority of these studies have focused on 
the evaluation of these novel forms of molecular nitrogen. However, Wang et a/.[3l| has 
described the possible synthesis route for all-nitrogen systems. Dixon et al. 32| has success- 
fully demonstrated the ability of the theory to predict the stability of such compounds. In 
spite of these efforts, so far no higher members (n>4) of polynitrogen family in free space 
have been synthesized experimentally . The present difficulties in synthesizing the polyni- 
trogen molecules motivates us to look for new avenues to realize the success of polynitrogen 
molecules as a possible HEDM. 

Further, the new polymeric forms of nitrogen have also been explored in extreme con- 
ditions of high pressure and temperature. The application of high pressure to nitrogen 
system leads to destruction of covalent bonds and may lead to formation of intermediate 
polymeric network of single and double bonded atoms which has been verified experimen- 



tally [33, 



34 



35, 



36|. Using ab initio calculations, other forms of polymeric nitro gen such 



as Black Phosphorous (B 
Cubic Gauche (CG) form 



'l(37|, A7 



, metallic 39], simple tetragonal phase [39j and 
36j have been proposed at extreme conditions. The CG form has 
also been observed experimentally 40j. Therefore, the success of ab initio studies 41] in 



predicting new polymeric phases of nitrogen in extreme conditions motivates us to study 
polynitrogen in other physical environment such as in confinement. 

More recently in 2008, nitrogen has been predicted to exist in polymeric form (iVg) in 
confinement of nanostructures i.e. Carbon Nanotubes. They proposed N 8 stabilizes as a 



chain which is stable up to room temperature [421 ] . Therefore, apart from extreme con- 
dition of high temperature and pressure, confinement offers an alternate environment for 
stabilization of polynitrogen com pou nds. Since the environment of Cqo has been extensively 



studied by our group [43, 



44 



45 



4a | as well as others 47], |48|] it offers a new possibility for 



polynitrogen encapsulation. The Cqq cage has already been explored for endohedral doping 



with alkali metals 



49j, transitional metals [50], non-metals 5lj and noble gases 52j. Among 



all endohedral complexes, nitrogen-doped C 60 are quite interesting and needs further inves- 
tigation. As per our knowledge, N@C m is the only endohedral complex which have been 



explored extensively theoretically 



47j, 



48 



53 



551 ] as well as experimentally [56(. However, 
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N 2 @Cqq was successfully synthesized using pressure heating j57| and one isolated study of 



its stabilization energy has been recently reported [58(. Till date, the higher members of 
nitrogen family (n>2) have not been explored at endohedral position in Cqq. New techniques 
have been developed to deal with nitrogen (as impurity) in nanostructures such as chemical 
vaporization deposition (CVD) method 60] which is capable of synthesizing nitrogen doped 
carbon nanotubes. And recently in 2008, C57N3 was reported using surface- catalyzed cy- 



clodehydrogenation from aromatic precursors 



61] and further expressing a strong possibility 



of synthesizing other nitrogen-based hetero and endohedral fullerenes. Keeping this in view, 
further efforts are needed to be invested in both experimental and theoretical sides to help in 
designing new ways for generating N n compounds and understand their stability in different 
environments. 

In this paper, we report the results of our ab initio calculations on N n @CQo for n<16, 
which describes the novel possibility of stabilizing the higher members of all-nitrogen 
molecules in confined environment of Cqo- We intend to show that 6*60 can act as an ideal 
candidate to trap nitrogen in polymeric form and can be used as a possible high energy den- 
sity material. The present study would provide some valuable information for synthesizing 
more stable all-nitrogen clusters to develop the novel materials for future. 

2. COMPUTATIONAL DETAILS AND VALIDATION 



We have used the Spanish 



(SIESTA) computationa 



tional theory method 



63 



64 



nitiative for Electronic Simulation with thousands of atoms 



code 62] which is based on numerical atomic orbital density func- 
651 ] . We have used Is 2 , 2s 2 , 2pl, 2py, 2p\ {AS^/z) configuration 



of atomic nitrogen corresponding to spin quartet state which is consistent with the experi- 
mental observation [56] . The calculations are carried out using generalized gradient approxi- 
mation (GGA) that implements Becke gradient exchange correlation functional (BLYP) 66]. 
Core electrons are replaced by non-local, norm- conserving pseudopotentials factorized in the 
Kleinman-Bylander form 67J , whereas valence electrons are described using linear combina- 



68 



but generalized 



tion of numerical pseudo atomic orbitals of the Sankey- Niklewski type 
for multiple-^ and polarization functions. In this work, we have used a split valence double-^ 
(DZP) basis set with energy shift equal to 100 meV. The real space cutoff grid energy is 
taken as 200 Ry. Periodic boundary conditions have been considered using simple cubic shell 
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of 2 nm that is large enough to avoid any significant spurious interactions with periodically 
repeated images. 

We first performed test calculations on the C 60 and N 2 molecules. In C 60 , we found the 
C=C and C-C bond distances as 0.140 nm and 0.146 nm in agreement with experimental 
values 69] . The ionization potential and electron affinity were found to be 6.9 eV and 2.70 
eV, in good agreement with experimental values of 7.5 ± 0.01 [3] and 2.689 ± 0.008 71] 
respectively. In N 2 molecule, the N-N bond distance and ionization potential are found to 
be 0.112 nm and 15.43 eV which are in good agreement with experimental values of 0.115 
nm and 15.60 eV respectively 721 ]. 

We now proceed with the investigation of nitrogen clusters inside Cqo. We started with 
different initial atomic configurations of N n for n = 1-16. The following considerations 
were taken into account while starting and determining the final structure for particular 
configuration, (a) Nitrogen atoms were placed sufficiently away from the walls of Cqq so 
that there is no likelihood to form any C-N bond, (b) Nitrogen atoms were allowed to 
form bonds among themselves and make polynitrogen compounds, (c) Each structure was 
allowed to relax till forces on each atom converges up to 0.001 eV/A. (d) The structure with 
minimum energy is considered as the final stable polymeric form of nitrogen, (e) All the 
optimized nitrogen clusters inside 6*60 were also relaxed in free space at OK. (f ) The structures 
of iV n @C6o were also relaxed at room temperature (300K) to confirm their stability (g) The 
spin polarized calculations were performed on final stable polymeric N u @Cqo to determine 
the spin states (h) Also harmonic vibrational frequency analysis on all optimized structures 
was performed to further ensure their ground state geometries. 



3. RESULTS AND DISCUSSION 



The optimized ground state structures of polynitrogen compounds inside 6*60 are shown 
in Figure 1-3. We started our calculation with single nitrogen, the nitrogen was placed at 
different positions inside C 60 and the SCF total energy isosurface is explored throughout 
region inside Cqq. We found the atomic nitrogen at center as stable structure with small 
zero point vibrations. Since the atomic nitrogen is highly reactive so if we place nitrogen 
at a distance of more than 0.150 nm from center of cage, it has tendency to form a weak 
covalent bond of the order of 0.154 nm with two adjacent carbon atoms of hexagon-hexagon 
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interface of Ceo inner surface [3]. When two nitrogen atoms are placed inside Ceo at different 
positions, they form N 2 molecule which is found to stabilize at the center of the C^o- The 
resultant structure is N 2 molecule with N=N bond distance of 0.112 nm which is same as 
in free space. 

In N 3 @Cqo, N 3 in linear form with symmetry is found to be stable at the center 
with N-N bond length of 0.119 nm. However, the triangular structure D^h is found to be 
another stable isomer inside Cqq. Interestingly, most of initial N3 configurations including 
the combination of N 2 and single N at distance results in an isosceles triangular structure 
with bondlengths of 0.133 nm and 0.158 nm and bond angles 53.5° and 73.1° respectively. 
The energy difference between both the structures is of the order of 0.9 eV. N± is found to 
exist as two dimers of N 2 with bond length of 0.111 nm and placed at a distance of 0.226 nm 
symmetrically across the center of cage. The arrangement of two N 2 molecules on either side 
of center of Cqq may be explained on the basis of extra stability of N 2 and its reluctance to 
form lower order bond. The structure with and D 2 h symmetry are found to be the other 
stable isomers. However, in free state N± has been predicted to be stable in C 2 h — T, 
and D 2 h symmetries respectively [3]. N& is found to be stable at center in perfect pentagon 
shape having D 5 h symmetry with each bond length of 0.132 nm and bond angle of 108.0 
°. We tried various initial configurations of N5 (including dimers and single atom). It was 
observed that all the configurations lead to the formation of pentagon structure which may 
be due to confined environment of Ceo- 

For n > 6, different combinations of N, N 2 , N3 and higher polynitrogens were considered 
as separate substitutional identities in Ceo- We found interesting pattern in the arrange- 
ment of nitrogen atoms in polynitrogen compounds, which are very stable and are placed 
symmetrically across the center. The other observed isomeric structures, though stable, are 
found distinctly less stable by an energy difference of ~ 0.5-2.0 eV. The minimum energy 
configuration structure of N e comes out to be a boat-shaped cyclic ring with bond lengths 
of the order of 0.128-0.135 nm and bond angles 108.8° and 122.5 °. In free state Nq is found 
to be stable in hexaazadewarbenzene structure with C 2v symmetry which is also one of the 



stable isomers in confinement 



75|. 



N 7 is found to exist as an assembly of a pentagon with bond lengths 0.128-0.132 nm 
and a dimer having bond length 0.111 nm as two separate units. The presence of N 5 as an 
independent identity is responsible for lowering the energy of Ni@Cqq indicating its extra 
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stability . However, one closed structure involving N$ and N2 do exist as a close energy 
isomer as shown in Figurel. Ng is found to exist in two planar rings having distinct bond 
lengths of 0.127, 0.138, 0.145 and 0.138 nm which are connected together with a weak N-N 
bond of 0.157 nm. N 9 exhibits a stable structure with N 5 and A" 4 planar rings connected 
together with a bond of 0.148 nm. The bond lengths for both the rings are in the range 
0.127-0.150 nm. iVio forms a structure with two planar symmetrical pentagon rings on either 
side of center of Cqo (D 5 h symmetry) connected to each other with a single bond of 0.150 
nm. Nu polynitrogen compound is also stable in ring structure with one planar pentagon 
and one non-planar hexagon connected by a bond distances 0.142 and 0.145 nm respectively. 
N12 stabilizes itself symmetrically into two non-planar boat shaped hexagons connected by 
N-N bond of 0.139 nm. 

For n>12, the polynitrogen compounds do not show any symmetry in distribution of 
nitrogen atoms due to the availability of less space. They make complex yet stable structures. 
iVi3 is found to exhibit a geometry with one planar pentagon connected to a non-planar closed 
N 8 ring with significant distortion on 6*60 cage. N14 and Ni 5 compounds make distorted 
endohedral complex with nitrogen undergoing chemisorption at the surface breaking the 
surface C-C and C=C bonds and making C-N bond of 0.136 nm. This leads to the formation 
of heptagon in the cage which results in very large distortion and the maximum and minimum 
diameter is observed to be 0.760 nm and 0.730 nm respectively. At n=16, the chemisorption 
of three nitrogen atoms results in the opening of the 6*60 cage. 

The encapsulation energy is defined as 

AE = E[N n @C m ] - E[C m ] - E[N n ], (1) 

where E[N n @C eo ] is the energy of the complex, i?[C 60 ] is the energy of an isolated C 60 and 
E[N n ] is the energy of an isolated N n cluster. These energies have been calculated for all 
the structures and are tabulated in Table I. We would like to add here that the energy of 
N n cluster in isolation is different from the N n in confinement as their structures are very 
different. These constrained polynitrogen molecules disintegrate into smaller clusters when 
relaxed in isolation. The encapsulation energy of N2@Cqo is found to be 13.8 kcal/mol, 



which is in agreement with the value of 9.11 kcal/mol reported using B3LYP [58|]. However 
we would like to mention here that B3LYP does not describe the non bonding systems 
accurately [59]. 



The encapsulated polynitrogen clusters results in distortion of icosahedral symmetry of 
the Cqo thereby altering the electronic properties of N u @Cqo significantly. The HOMO- 
LUMO gap variation as a function of number of nitrogen atoms shows gradual decrease 
with higher values for n= 2, 4, 6, 8, 10. To visualize this we have plotted the Kohn Sham 
energy levels for N u @Cqq complexes using DZP basis in Figure 4. Also electronic density 
of states (DOS) for the structures are shown in Figure 5-6. We carried out spin polarized 
calculations on all the optimized geometries. Multiplicities defined as (2s+l), where s is the 
spin of the system, are calculated from the difference of spin up and spin down electrons (Q 
up down) and are tabulated in Table II. 

For N@Cqq we found a spin quartet state when the N atom is placed at the center of 
the Cqq cage i.e N atom retains its atomic nature, which is consistent with the experimental 
ESR or EPR spectra of N@Cqq 47J, |53|, |54| • However when the N atom moves toward the 



surface and makes a covalent bond with the cage, it shows a spin doublet state consistent 

n 

with experimental results [73J . All the structures with even number of nitrogen atoms show 
spin doublet states (s=l/2) whereas when odd number of nitrogen atoms are present they 
show spin singlet (s=0) state except for n=13 which shows a spin quartet state (s=3/2). 

The net charge on the polynitrogen compounds and the fullerene cage atoms are investi- 
gated by Mulliken charge analysis. We would like to mention here that charge transfers in 
SIESTA are dependent on the choice of basis set. Therefore, we carried out test calculations 
using SZ (single zeta), DZ (double zeta) and DZP (double zeta with polarization) (with 
energy shift ranging from 50 meV -250 meV) basis sets on N@Cqo- It was observed that 
using SZ basis N shows a small charge loss of 0.04 electrons to Cqq cage, whereas using DZ 
and DZP basis N shows a small charge gain from the cage. The magnitude of charge transfer 
is found to be higher when DZ is used compared to that using DZP. Within DZP basis set 
as we increase the energy shift the mulliken charge on N tends to increase marginally from 
0.027 (energy shift 50 meV) to 0.040(250 meV). Our work on azafullerenes using DZP basis 
set produced reliable charge transfers in C6o- n N n [ifl and it reproduced the experimental 
parameters of Cqo and N2 very accurately. 

For N@Cqo, the small charge transfer of 0.04 electrons from the cage to N at center is 
found to be consistent with spin quartet state as obtained using spin polarized calculations. 
However, there are conflicting claims regarding charge transfers in N@Cqq using various ab 
initio calculations. Kobayashi et al. Tj\ using 6-31G(d) basis set and B3LYP functional has 



S 



reported slightly large spin density at N in N@Cqq than in iV@CVo consistent with observed 



hyperfine coupling constants. Lu et al. 



781 ] using user defined basis set (considering 2s and 



2p for C and N) and carrying discrete variational local density functional calculations had 
reported charge transfer of 0.11 from C^q cage to N atom at the center. Further. Mauser et 
al. j3] using semi empirical and DFT(B3LYP) and Greer [79J using TZVP spin-unrestricted 
DFT optimized geometry within DZP basis sets have predicted no charge transfers between 
N and 6*60- For n=2-10, we found charge transfer from cage to N n , which may be explained 
on the basis of higher values of electronegativity and ionization potential of N than C atoms. 
The Mulliken charge calculation suggests electron charge transfer ranging between 0.05 to 
0.31 from carbon atoms to the polynitrogen compounds inside up to n=10 except for N 2 
which shows inverse nature. For n>10, charge transfer of 0.05-0.47 electrons takes place 
from nitrogen to carbon atoms of fullerene. This reversal can be attributed to availability 
of less space and overlapping of orbitals which may be further responsible for destabilizing 
the polynitrogen compounds beyond n=13. As per our information, there are no references 
for comparison to any theoretical work on polynitrogen clusters in Cqq till date. 

From the structural analysis of all the N n @CQo structures and their isomers it may be 
concluded that N atoms have preferential sites for its stabilization inside fullerenes. More- 
over, due to confinement the nitrogen atoms are forced to form lower order bonds. We found 
that the threshold size of polynitrogen compounds to exist inside as stable molecule is A^ 13 . 
Any further addition of nitrogen inside leads to the formation of C-N bond and heptagon 
involving nitrogen atom on the surface. All the nitrogen clusters in N h @Cqo for n>4 has 
70-80 % single bond character which confirms the polymeric nature of nitrogen compounds 
as shown in Figure 7. The structural distortion expressed in terms of elongation defined as 
e = (-^r- — 1) where R and Rq are the average radii of expanded and pure Cqo is shown in 
Figure 8. The percentage elongation shows a linear increase with the number of nitrogen 
atoms up to except for N2 which shows small contraction in 6*60 cage. The maximum 
threshold elongation is found to be ~ 6% for iV 15 . Indeed, this elongation leads to estima- 
tion of extra pressure. Successive addition of nitrogen atoms exert this extra pressure from 
within when their number is increased beyond 2, till the stability limit of C^o is reached. 

Further, the stability of each structure is verified by displacing the nitrogen atoms up to 
0.5 nm and observing the result. We find that after relaxing the structure , it regains its 
optimized form. Moreover, harmonic vibrational frequency analysis is carried out on all the 
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ground state structures and the absence of any imaginary frequency value further verifies 
the true energy minima. The ground state geometries for polynitrogen N n compounds, 
when relaxed in free space, leads to fragmentation with large number of N 2 molecule units. 
This suggests that the polynitrogen compounds are forced to remain stable inside due to 
confinement and any mechanism triggering the breaking of Cqo cage may lead to the release 
of enormous amount of energy. The energy released can be approximately estimated based 
on the energy of N n Ceo and 6*60 and equivalent (|n of N 2 molecules). Assuming all n atoms 
of N in N n Ceo clusters are singly bonded (with energy 158.9 kJ/mole), the energy released 
on formation of N 2 molecules comes out to be around 3-6 kcal/gm of nitrogen atoms. This 
is several times larger than the energy release in the case of TNT. 

We performed constant-temperature constant-volume molecular dynamics for iV n @C 6 o 
at a temperature of 300K with time period 10 ps with interval of 1 fs. It was observed 
that for n <12 the polynitrogen structures remain intact and do not show any significant 
deviation from their structures at OK which demonstrates the stability of polymeric nitrogen 
structures. But for n >13, the 6*60 cage opens up due to formation of C-N bonds as a result 
of thermal vibrations. Interestingly, our preliminary calculations on addition of few water 
molecules in iV n @C6o for n=8-10, indicates that the polynitrogen dissociates resulting in 
breaking of 6*60 cage. 

The experimental synthesis of these novel hybrid systems would be a challenging task. 
Owing to the instability of N n clusters in free space, the most feasible and promising way 
is to trap them in some nanostructures such as fullerenes, nanotubes and transition metal 
complexes. It appears to us that N u @Cqo could be synthesized by nitrogen enrichment using 



prolonged ion implantation of N ions or pressure heating (57|. It is well known that laser 
ablation of graphite can produce copious amounts of Cqq. When graphite, intercalated with 
desired specie of atom/ molecule is laser ablated, endohedral Cqq results. Therefore, even if 
short-lived N n clusters could be intercalated before laser ablation, there is good possibility 
of their getting trapped inside Cqo cage. 

4. CONCLUSIONS 

The polynitrogen compounds formed inside Cqq are found to be quite stable in confine- 
ment. The energy difference of the optimized ground state geometries to the isomeric states 
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is found to be significant, thereby indicating clear choice of ground state structures. The 
spin states of N n @C eo are found to be quartet for n=l and n=13, doublet for other odd 
n and singlet for even n. The reversal of polarization and the space available inside C 60 
may be responsible for determining the amount of nitrogen encapsulation inside fullerene. 
The maximum number of nitrogen atoms that can be encapsulated inside Cqq is 13, which 
can form stable structure at OK. The endohedral molecules N n @Ceo for n < 12, are found 
to retain their structure at room temperature. This suggests the strong possibility of 
synthesizing polynitrogen compounds in confinement which could be triggered with suitable 
combinations of oxygen and hydrogen for extracting tremendous amount of energy. The 
amount of energy density from conversion of lower order bonds to triple bonds comes out 
to be about three times in contrast to the conventional chemical energetic materials. 
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Table I: The encapsulation energy in eV defined as delta E = E comp i ex -Ec 60 -E^ n , 
where £V n is the energy of the N n compounds in isolation in free space. Note the 
energy considered here are of theoretically lowest energy structures which are very 
different from structures in confinement and do not exit in free space for n>3. 



N n <i 


5)C 60 


Total Energy (eV) 


Encapsulation energy (eV) 


NCQ 




-9490.2 


-2.4 






-9766.0 


-0.6 






-10030.5 


0.2 






-10301.0 


-3.5 


N 5 1 




-10566.7 


2.9 






-10831.2 


11.8 


N 7 1 




-11097.8 


10.7 


N 8 Ci 




-11360.1 


9.8 


N 9 Ci 




-11624.6 


16.4 


N l0 ( 




-11889.4 


25.5 


N n ( 


ac 60 


-12152.6 


33.9 


N 12 ( 


ac 60 


-12414.9 


37.7 


N 13 ( 


ac 60 


-12678.6 


39.2 



17 



Table II: The Spin multiplicites of N n @Cm complex and difference between spin up 



and spin down charge densities 
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(a) 0>) 




(e) (f) 



Figure 1. Geometries of N n @Ceo (a) N2 placed at the centre of C60 (b) N3 as a linear chain placed 
symmetrically (c) N4 exists as two independent Ni molecules placed symmetrically across centre (d) & 
(e) N5 and N6 exist as a pentagon ring and a non-planar distorted hexagon with 100% single bond 
character (f) N7 exists as N5 and N2 as separate units, however, its nearly isoenergetic structure is shown 
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(a) (b) 




(e) (f) 

Figure2. Geometries of N n @Cso (a) Ng exist as two planar rings of N4 placed symmetrically and 
connected by weak N-N bond (1.57 A) (b) N9 exist as two planar rings of N4 and N5 connected by weak 
N-N bond (1.48 A) (c) N10 exists as two planar and slightly inclined pentagon rings connected by weak 
N-N bond (1.50 A) (d) Nn exist as interconnected structure of planar N5 and non-planar N& placed 
across centre (e) N12 exist as symmetrically placed two boat shaped hexagons connected by N-N bond 
distance of 1.39 A (f) No exists as interconnected structure of planar N5 and a heavily distorted Ng 
however, the less interstitial space inside Cf,o forces the nitrogen atoms to make bonds with surface 
carbon atoms. The above structures of polynitrogen molecules in C60 show nearly 70-80% single bond 
character. 
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(c) 



Figure3 (a) Ni4@C60 exits as a distorted complex with N12 existing as interconnected structure of N5 and 
N7. Two nitrogen atoms undergo chemisorption resulting in 62 atom cage with the formation of two 
heptagons on surface, (b) Ni5@C60 doesn't show any pattern in structure of polynitrogen molecule. 
Hexagons and pentagons of Ceo show heavily distortion with C-C bond distances increasing up to 1.65- 
1.67A which is the breaking point for C60 cage, (c) Ni6@C6o exists as a polymeric complex with three 
nitrogen atoms getting substituted at the surface and making a heptagon and octagon which henceforth 
leads to the opening of the cage indicating the threshold limit for nitrogen encapsulation at OK. 
In above structures, due to less available space nitrogen tends to make bonds with the surface atoms as 
well as undergoing chemisorption. These structures are unstable at room temperature. 
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Figure 4. Shows the Kohn-Sham energy eigen values for N 
clusters (n=1 ,2,5,6,8,9,1 0,1 2,1 3) 
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